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I review the progress made in the calculation of the electroweak penguins Q7,s in the framework 
of the Hadronic Approximation to large- N c QCD, and give an update of results. 



A nonzero value of e'/e signals the existence of direct CP violation in the decay K — > tttt 
and is a measure of the amount of CP violation in the Standard Model. At present the world 
average is given by a 



- = (16.6 ± 1.6) x 10" 



(1) 



To compute this number in the SM has become a major theoretical challenge!. Roug hlyi, e'A 
is dominated by the difference of K — > tttt matrix elements of the Qq and Qs operators, 



^{6,8} = - i 8 ' 12 M ( S LQR) ((lR d L) , 



(2) 



q=u,< 



where the factors e q = (2e/3, — e/3, — e/3) are the quark electric charges. These factors make a 
profound difference between Qq and Qs at the meson level as they change their chiral properties 
making Qs & n operator which starts at 0(jP), whereas Qe only starts at 0(p 2 ). This makes Qs 
a simpler object to deal with and allows a connection to the well-known problem of the tt + — tt° 
electromagnetic mass difference. 

Strong interactions make these operators change with the renormalization scale and mix 
with others so that, strictly speaking, they can never be considered in isolation. In particular 
Q§ mixes with Q7, 

Q 7 = 6 e <? (slJ^cIl) {q R -i^q R ) ■ (3) 

q=u,d,s 

This apparent complication will actually be very helpful in the calculation of Qs- 



The main effect is due to the fact that Q7 and Qs contribute to the operator 



Tr(u\f>rfQ R ) , (4) 

where (A^ )y = ^3^2 and Qr = diag(2/3, — 1/3, — 1/3). This operator is a "flavor-rotated" 
version of the operator responsible for the ir + — 7r° electromagnetic mass difference, which is 
a problem that has been reasonably understood for a long time El. This relationship is the 
underlying reason why both the physics of the 7r + — 7r° mass difference and the bosonization of 
Qts into the operator (||) are governed by one and the same Green's function, which turns out 
to be the analog of the vacuum polarization but between the left- and right-handed currents 
Ql,r1^Ql,r- Furthermore, Qj can be looked upon as due to the exchange of an imaginary 
"photon" with couplings to the quarks which are nondiagpnal in flavor and whose "propagator" 
is unity instead of the usual 1/q 2 . One then obtains that Q 

Q 7 = 6 < OiOu) > Tr (uX^U^Qr^ , (5) 

with 

< Oi(/i) >=< s L ^d L dR^SR >= ^ (^-Mg^ J ^^n^(g)^ , (6) 

where current conservation implies il^(g) = {(f , q v — q 2 g^ u )H-LR(l 2 ), since we are in the chiral 
limit. Because our imaginary photon has unity as its propagator, this integral is divergent and 
has to be regularized and renormalized. The subscript MS refers to this fact. Eq. (|6|) shows 
an example of how a coupling constant in the effective meson Lagrangian (i.e. < Oi(/i) >) is 
related to integrals over euclidean momentum of QCD Green's functions. This is also true in 
general. 

However, in order to really compute < 0\ > one still needs to know the function HLR(q 2 ) for 
all values of q 2 which is not known. The large- N c limit of QCD simplifies the problem because 
it determines the analytic structure of this function: it has to be meromorphic. Therefore it 
can only have poles (and no cut), which in physical terms correspond to the meson states. 
However this is still not enough as the number of poles is infinite and, since the solution to 
large- N c QCD has not been found, the location of these poles and their residues are unknown. 
It is at this point that a further approximation has to be made. This approximation, which 
we name "The Hadronic Approximation" to large- N c QCD, is defined by keeping only a finite 
number of resonances, whose residues and masses are fixed by matching to the first few terms 
of both the chiral and the OPE expansions f| of Hlr(Q 2 ) fl- This is called in mathematics 
a rational approximant and is an interpolating function, ratio of two polynomials, which by 
construction has the same low- and high-Q 2 behavior as the full I^lr{Q 2 ) of large- QCD. This 
approximation rationalizes and systematizes old phenomenological approaches such as Vector 
Meson Dominance; incorporating chiral symmetry (to control long distances) and the OPE 
(to control short distances). It is also an improvement over approaches where only the OPE 
constraints were considered, such as QCD sum rules, because it also uses chiral symmetry at 
long distances. 

An analysis of Aleph datai shows that just the pion and one vector and one axial- vector 
states do a pretty good job provided their masses and decay constants fulfill the above-mentioned 
long and short-distance constraints. I urge the reader to take a look at the curves in Fig. 1 of 

a In the past we used the term "Minimal Hadronic Approximation" for cases in which only the leading term in 
the OPE expansion was used in the matching. 

^Notice that this is possible because Hlr(Q 2 ) has no parton-model log Q 2 ; it would not be possible in the case 
of the vector- vector correlator, for instance. 



the first paper in Ref. 0. In this case one obtains the remarkably simple expression 



- Q 2 Klr{Q 2 ) - (Q2 + M^ + Mi) ' (7) 

where f = 87 ± 3 MeV, My = 748 ± 29 MeV and M v /M\ ~ 0.50 ± 0.06 1. These values do 
not have to be exactly equal to the physical ones although it is natural to expect that, if the 
Hadronic Approximation and large- N c work, they have to be close to each other. 

Having the function Ulr(Q 2 ), one can plug it in Eq. (||) and compute the integral. The 
result is 



<Ox>(M) = -*L log 



32vr 2 ! _ K 



My 

M\ i M 2 



f 2 



(8) 



with A 2 = /u 2 exp(l/3 + k), where k = —1/2(3/2) in NDR(HV) schemes, respectively. The 
control over the renormalization scheme (75, evanescent operators, etc..) is the result of the 
OPE constraints used in the construction of the Hadronic Approximation □. It is obvious that 
knowing < 0\ > is tantamount to knowing any matrix element like, e.g., < tttt \ Qj \ K >: one 
only expands the U in Eq. (||). 

The calculation of the matrix elements of Qs is a bit more tricky because it is the product of 
two scalar-pseudoscalar densities -see Eq. (||)-. The previous analogy with an imaginary "pho- 
ton"! still works and leads to a result similar to (||),(§D but involving these scalar-pseudoscalar 
densities, to witi 



Q 8 = -12 < 2 (p) > Tr ([/Af%tQ fl ) T , (9) 



where 



< 2 (» > = < s L s R d R d L > 




~<^Hb+| /ttttTT /' /J ''- W[>1l<U<U-) s«*lW |0)r } (1.0) 



MS 



Taking the large- N c limit selects the quark condensate in the previous equation. However, there 
is some evidence -although circumstantial B- that this limit in the scalar-pseudoscalar sectox 
may have somewhat large subleading corrections. This is why in our original paper in Ref. Q 
we decided to take a short detour. The operator Qj mixes through the renormalization group 
into Qg; this is tantamount to saying that < 2 > controls the large-Q 2 fall-off of Hlr(Q 2 ). 
Consequently, by expanding Eq. (^) and selecting the coefficient in 1/Q 6 , one obtains 

-Q 6 n LR (Q 2 ) ^ flM v M\ (fromEq.©) 



OL 



—Q Ulr{Q ) q tt^ 16to s \1 + < 2 (ji) > + ■ ■ ■ (from the OPE) , (11) 

where £ = (25/8,21/8) in the (NDR, HV) schemes !0, and the ellipses stand for numeri- 
cally negligible terms ^. Inputting the value of a s (2 GeV) ~ 0.33 ± 0.04 one easily obtains 
< 2 (2 GeV) > from matching the two asymptotic behaviors in Eq. (|TT|) . Saturation of this 
value for < 2 (2 GeV) > by the quark condensate in Eq. ( pTj| ) leads to {ipip)j]-g(2 GeV) ~ 
-(300 ± 20 MeV) 3 , where the error is a combined estimate of unknown 1/N C and a s correc- 
tions. 

At any rate, having the value of < 2 (2 GeV) > one can easily compute any K — ir matrix 
element through Eq. (^). Table 1 shows our results together with those obtained by other 

c This "photon" is now even more estrange: it doesn't have spin one!. 

d At the time we wrote our paper □ the term proprtional to £ was not known. 



Refs. 


M 7 (NDR) 


M 7 (HV) 


Mg (NDR) 


M 8 (HV) 


AfThis worki Knarht pt nM 
NarisorO r, 
Cirigliano et-alu 
Maltmanlii, 
Bijnens et a^Eji— i 
Battacharya e£i-«/ll3 
Donini et alfJ 


nil ± o 03 
0.17 ±0.05 
0.16 ±0.10 
0.21 ±0.03 
0.24 ±0.03 
0.32 ±0.06 
0.11 ±0.04 


67 ± 90 

0.49 ±0.07 
0.46 ±0.08 
0.37 ±0.08 

0.18 ±0.06 


^.J4; _1_ \J. 1 O 

1.4 ±0.3 
2.22 ± 0.67 
1.65 ±0.45 
1.2 ±0.9 
1.2 ±0.2 
0.51 ±0.10 


Q W + (17Q 

Zi.iJZj _L_ Lf. 1 J? 

2.46 ±0.70 
1.84 ±0.46 
1.3 ±0.9 

0.62 ±0.12 



Table 1: Summary of results for M 7 ,8 s< (7r7r)/ =2 |Q7,8|A' > (2 GeV), in units of GeV 3 . 



groups. A look at the central values (and errors!) in Table 1 reveals that, even though a lot of 
progress has been made, meson weak matrix elements are a very hard problem to understand 
and we are not quite there, yet. Although it should be obvious, it is worth remarking that lattice 
and continuum approaches may nicely complement each other in this effort. 

I would like to finish by emphasizing that the Hadronic Approximation to large- N c QCD is 
not limited to simple operators like Qt^ only, but isa general framework which has also been 
successfully used in a variety of other problems: -Ba'"> ^° — ► e + e~ and 77 — > /x + /i~ll3, hadronic 
contributions to g — 20, etc.; even in those cases when the underlying QCD Green's function 
could not be related to any experimental data -which are the most-. More work is in progress. 
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